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Introduction. 


vjr The  main  objective  of  the  present  note  is  to  solve  a  problem  proposed 
by  Weide  (1978)  concerning  the  complete  convergence  of  certain  random 
variables  associated  with  Karp's  probabilistic  analysis  of  the  traveling 


salesman  problem  (Karp  (1976),  (1977)).  ^ 

To  set  the  problem  precisely  let  X^,  1  <  i  <  «,  be  independent 

2 

random  variables  uniformly  distributed  on  the  unit  square  [0,1]  ,  and 


let  Tn  denote  the  length  of  the  shortest  path  (in  the  usual  Euclidean 
distance)  irtiich  connects  each  element  of  . . .  ,Xh). 

It  was  proved  by  Beardwood,  Hal ton,  and  Hammersley  (1959)  that 


lim  Tn//n  =  P 

n  oo 

with  probability  one  for  a  finite  constant  P.  This  fact  was  central  to 
the  motivation  behind  Karp's  algorithm,  but  as  Weide  (1978)  points  out 
the  Karp  algorithm  actually  calls  for  the  following  stronger  result  to 
be  proved  here: 

Theorem  1.  There  is  a  constant  P  such  that  for  all  e  X),  one  has 

on 

E  P(|T  //n-P  |  >  c)  <  «  . 

n=l 

This  type  of  convergence  is  usually  called  complete  convergence,  and 
Theorem  1  stands  in  a  similar  relation  to  the  Beardwood,  Hal  ton,  Hammersley 
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Theorem  as  the  Hsu-Robbins  Theorem  stands  in  relation  to  the  strong  law 
of  large  numbers  (Lukacs  (1968),  Hsu  and  Robbins  (19^7)).  The  "easy -half" 
of  the  Bor el -Cantelli  lemma  shows  that  Theorem  1  implies  the  Beardwood- 
Halton-Hammersley  Theorem  and  the  "hard-half"  of  the  Borel -Cantelli  lemma 
shows  how  Theorem  1  is  necessary  in  modeling  contexts  where  problems  of 
increased  size  are  generated  independently  of  previous  problems.  (For  a 
full  discussion  of  independent  versus  incrementing  models  for  random 
problems  one  should  consult  Weide  (1578)). 

The  proof  of  Theorem  1  is  given  in  the  next  section  and  depends  upon 
shapening  a  subadditivity  argument  which  has  been  useful  in  more  general 
contexts  (Steele  (1979))*  The  third  section  discusses  a  generalization 
of  Theorem  1. 


2.  Proof  of  Theorem  1. 

Let  N  =  N(t)  denote  a  Poisson  counting  process  with  constant  growth 

p 

rate  1.  Also,  for  any  Ac[o,l]  let  T^(A)  denote  the  length  of  the 
shortest  path  through  the  points  A  n  The  method  of 

proof  rest  upon  developing  recursions  for  the  functions 

l/2 

<p(t)  =  E  Tt([0,l]2),  *(t)  =  (E  T2([0,1]2))  ,  and  V(t)  =  *(t)-<p2(t). 

First  one  notes  the  following: 

Lemma  1.  There  is  a  constant  c  such  that  for  any  set  i  c  [0,1]  with 
n  elements  there  is  a  path  of  length  no  greater  than  c/n  through  i. 

This  lenmta  is  easeily  proved,  but  for  a  proof  which  yields  a  good 
value  of  c  one  cam  consult  Few  (1955). 
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Now  let 


2 

s  =  Tt([0,|]  )  +  Tt([|>l]2)+Tt([|>l]x  [0,|])  +  Tt([|>l]X  [0,|]) 
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and  note  by  elementary  geometry  that  ([0,1]  )  <  S  +  4.  By  well-known 
properties  of  the  planar  Poisson  process  one  has  that  the  four  summands 
of  S  are  independent  and  identically  distributed.  By  scaling  one 
also  notes  E  Tt ([0,|]2)  =  |  <p(t/4 )  and  (E  T^([0,|]2)  )^2  =  |  *(t/4). 
Thus,  taking  expectations  of  (S+4)  one  has 


(2.1) 


\|r(t)  <  f  (t/4)  +3  qT(t/4)  +l6<p(t/4)  +  16  . 


To  simplify  (2.1)  let  ^(t)  =  \|r(t2),  <p^(t)  =  <p(t2),  and  V^(t)  =  V(t2) 
and  note  that  (2.1)  implies 


*^(2t)  <  *2(t)  +  3  (t )  +l6^(t)  +  16 


which  entails 


V1(2t)  <  V1(t)  +  4  <p2(t)  -  qpj (2t )  +16  ^(t)  +  16 


and  dividing  by  (2t)2, 


(2.2)  V1(2t)/(2tf  -  (lAJV^tJ/t^  T2(t)/t2  -  q>2(2t)/(2t)*  +  U (q^  (t  )+l )/t; 


2  M-l 

Applying  (2.2)  successively  to  t,2t,2  t,...,2  t  and  summing  yields 


£  V  (2kt)/(2kt)2  -  J-  Y  V  (2kt)/(2kt)2 
k=l  k=0  A 


<  <f>2  (t  )/t2  -  <p2  (2Mt )/  (2Mt  )2  +  4  (<p  (t)+l )  (  £  2“2k)t‘2 
~  -1-  A  x  k=0 


and  consequently 


(2.3)  I  V1(2kt)/(2kt)2  <  5’(V1(t)/t2+«p^(t)/t2+8(<P1(t)+l)/t2)  <  » 


In  terms  of  V(t),  (2.3)  becomes  for  u  =  tc 


(2.4) 


uu  ,  , 

£  V(4  U.)/ (4  u)  <  00 
k=0 


Now  by  Beardwood-Halton-Hammersley  and  dominated  convergence  theorems  we  have 
cp(4ku)/  (4ku)1//2  as  k  -*oo.  Hence  by  (2.4)  and  Chebychev’ s  inequality, 

(2.5)  £  P(|T  v  ([0,l]2)/(4ku)l/2-P|  >€)<». 

k=0  4  u 

Now  for  all  n  >  4P  there  is  a  4P  ^  >  m  >  4P  such  that  4  m  <  n  <  4  (m+1 ) 
for  some  k  >  0.  Hence  by  Boole's  inequality  and  the  mo no tonicity  of 

Tt[0,l] 

(2.6)  £  P(|T n([0>l]2)/n1/2-e|  >  e) 
n  0 

l^p+1 

<1  L  p(T  v  >  (P+e)(4km)l/2) 

“k=1m^  * 

4p+l 

♦  I  £  P(T  k  <  0-e)(4k+1(»+l)l/2) 
k=1  m=4p  14  m 
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»  4P+1  ,  -1/2 

<£  l  p(Tk  /(*♦  (»+D  '  >  (P+e)(l+4-p)  ) 

k=l  4  (m+1) 

~  **P+1  k  1/2  n  1/2 

+  I  I  P(T  k  /  (k  ®)  <  0-e)(l+4-p)  )  . 

k=1  m=4p  U  “ 

-D  -1/2 

3y  choosing  p,  one  can  guarantee  that  (1+4  y)  (P+e)  >  P  and 
(l+r_^)(P-e)  <  P  so  (2.5)  insures  that  the  last  two  sums  in  (2.6) 
converge,  so 


(2.7)  £  P(|T  ([0,l]2)/nl/2-p|  >  o)  <  oo  . 

n=0  n 

To  obtain  Theorem  1  from  (2.7)  we  note  that 

(2.8)  P(|T  «0,l]2)/nl/2-  P(  >  c)  *  £  P(|Tm/nl/2-0|  >  e) 

m=0 

oo  ,  /_  m  -hi 

>  I  P<Vn  /  £  s*‘>  •  Tf- 

m=n 

n  .  z*  m  -m 

+  £  P(T  /nl/2  <  P-e)  •  2«_ 
m=0  “  mJ 

>  P(Tn/n1/2  >  P+€)(J«)  +  P(Tn/nl/2  <  P-e)(J-)  . 


In  the  last  inequality  above  one  uses  the  fact  that  P(Tr  >  \ )  is 

monotone  increasing,  and  the  fact  that  £  nme”n/ml  and 

m  >  n 

£  nme“n/®t  each  exceed  r-  .  The  bounds  (2.7)  and  (2.8)  complete 
m  <  n 


the  proof  of  Theorem  1. 
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3.  Further  Results 


In  the  previous  section  the  aim  was  to  give  the  most  direct  proof 
possible  of  the  result  conjectured  by  Weide  (1978).  The  complete  conver¬ 
gence  proved  there  can  easily  be  extended  to  the  broader  context  of 
subadditive  Euclidean  functionals.  Since  the  method  of  proving  complete 
convergence  in  this  context  is  directly  parallel  to  the  preceeding 
argument,  it  is  sufficient  to  state  the  more  general  result. 

To  introduce  Euclidean  functionals,  let  L  denote  a  real  value 
function  defined  on  the  finite  subsets  |Rd,  d  >  2.  It  will  be  assumed 
that  L(X^,Xg, . . .  ,Xfl)  is  measurable  whenever  X^,  1  <  i  <  n,  are 
measurable,  and  that  the  following  four  assumptions  are  satisfied: 

AI«  L^x^^Xg,. . .  ,^x^ )  =  Q:l(x^,X2,.  •  *  ,x^)  for  Q-  ^  0  # 

A2.  L(x1+t,x2+t, . . . ,xn+t )  =  L(x1,Xg,. . . ,xr)  for  t  €  Bd  . 

A3.  L(Xl,x2,...,xn,xn+1)  >L(xlfx2,...,xn)  . 

A4.  Var(L(X1,X2^*  •  •  >Xn))  <  oo  when  X^  are  independent  and 

uniform  on  the  unit  cube  [0,l]d  . 

lere  is  one  further  assumption  -which  is  needed  and  it  is  the  only 
one  ldiich  is  not  trivial  to  verify  in  most  applications. 

Suppose  that  (Q^:  1  <  i  <  md)  is  the  partition  of  [0,l]d  into 
sub-cubes  of  edge  l/m  and  let  tQ^  =  (x:  x=ty,  y  €  Q^).  The  subadditivity 
assumption  is  the  following: 

A5*  Biere  is  a  constant  C  >  0  such  that  for  all  positive  integers 
m  and  positive  reals  t  one  has 

d 

d  m  d-1 

k({xl,x2>  •  •  *  A  [0,t]  )  <  £  * #  *  #^n^  ^  tQ^ )  +  Ctm  • 
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It  is  proved  in  Steele  (1979)  that  if  L  satisfies  the  preceding 
five  assumptions  and  (X^)  are  independent  with  a  uniform  distribution 
in  [0,l]d  then 

lim  L(X1,X2,...,Xn)/n(d_l)/d  =  0 

n  oo 

with  probability  one  for  some  constant  0  <  P  <  oo  . 

One  can  easily  check  that  the  preceeding  result  implies  the 
Beardwood,  Halton,  Hammersley  theorem.  The  main  observation  to  be 
made  here  is  that  this  result  also  can  be  sharpened  as  in  Section  2 
to  obtain  the  following: 

Theorem  2.  Suppose  L  is  a  functional  which  satisfies  assumptions 
A1-A5-  If  1  <  i  <  oc]  are  independent  and  uniformly  distributed 

in  [0,1]^,  then  there  is  a  constant  0  <  p  <  oo  so  that  for  all 

€  >  0 

Z  )/n^d"1^d-e|  >  e)  <  oo  . 

n^l 

There  are  many  functionals  which  arise  in  the  theory  of  algorithms 
and  which  satisfy  A1-A5.  The  easiest  examples  are  those  associated 
with  the  Steiner  tree  problem  and  the  rectilinear  Steiner  tree  problem, 
but  in  essentially  any  geometric  problem  which  deals  with  minimized 
lengths  one  can  find  natural  functionals  which  meet  the  conditions  of 
Theorem  2. 


7 


As  a  final  point  one  should  note  that  both  Beardwood,  Halton,  and 
Hammersley  (1959)  as  well  as  Steele  (1979)  contain  results  valid  for 
random  variables  with  non-uniform  distribution.  The  approximation 
processes  used  in  these  papers  to  extend  the  uniform  case  can  again  be 
applied  here  although  to  do  so  would  require  considerable  space.  Since 
the  algorithmic  applications  and  the  program  begun  by  B.  W.  Weide  (1978) 
are  ably  served  by  Theorems  l  and  2,  these  last  extentions  have  been 
omitted. 
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Let  X^,  1  £  i  <  00  ,  be  uniformly  distributed  in  [0,1]  and 
let  be  the  length  of  the  shortest  path  connecting  . 9X 

It  is  proved  that  there  is  a  constant  0  <  6  <  00  such  that  for  all 

e  >  0 

oo 

l  P(|T  /Vn-B|  >  e)  <  »  . 
n-1  n 

This  answers  a  question  posed  by  B.  W.  Weide. 


